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Abstract. Using a quantum electrodynamical approach, we derive the scattering 
phase matrices for polarized radiation in forbidden line transitions and in the presence 
O I of an external magnetic fields. The case of ( J = ^ 2 — >• 0) scattering is considered as 

an example. The non-magnetic Rayleigh scattering phase matrix is also presented. The 
C3 '. Stokes profiles in a single scattering event are computed for the strong field (Zeeman) 

and weak field (Hanle) limits, covering also the regime of intermediate field strengths 
(Hanle-Zeeman) . 
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cn ' 1. Introduction 

in 

^ . In astronomical spectroscopy, the so-called forbidden transitions are as important as the 
O allowed transitions p. They are seen [II S 11 El El E El l M [III M [El 112 [ISl llEl [I^ 

in the spectra of rarefied gases and plasmas under special conditions, such as are found 
^ . in some nebulae, solar corona, parts of active galactic nuclei and the extreme upper 
atmosphere of the earth. Historically, the term "forbidden" was associated with all 
those atomic transitions which do not involve emission or absorption of the electric 
dipole radiation. The forbidden lines may some times account for 90% or more of the 
total visual brightness of an object like a planetary nebula, so much so that, when 
such lines were first seen in the 1860's, they were thought to be due to a new element 
which was called as Nebulium. Forbidden lines are often unaffected by departures 
from local thermodynamic equilibrium (LTE) and are also insensitive to atmospheric 
temperature uncertainties as well as micro-turbulence. As such they have been quite 
useful to measure abundances, because at low densities of the order of a few atoms per 
cubic centimeter, the strength of a spectral line grows as A^^ where denotes the total 
number of atoms in the concerned atomic state, whereas at high densities it grows as A^. 
Scattering of radiation by atoms involves absorption followed by spontaneous emission 
and the scattering phase function has to be generalized [H] to a 4 x 4 phase matrix to 
include polarization state of the incoming or outgoing radiation. 
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With technological advances in spectropolarimetry, accuracies of the order 10~^ 
have been achieved in measuring the polarization expressed through Stokes parameters. 
Using such instruments, it has been possible to discover a wealth of new and unexpected 
phenomena that are taking place in the outer layers of the solar atmosphere [19]. These 
recent observations have opened a new window to the diagnosis of the polarized light 
scattering phenomena in the Sun and the stars [20l[ni[23[23l[2l[25l[26l[27l[28l[29]. It 
is well-known that magnetic fields are present almost everywhere in the universe |30j . 
Therefore we focus our attention on scattering in the presence of magnetic fields. In 
our previous papers we developed a quantum electrodynamical (QED) approach to 
line scattering in the presence of external electric and magnetic fields of arbitrary 
strengths [31], [32]. That approach can take account of all multipole type transitions, 
apart from the dominant one, namely, the dipole scattering. In this paper, we restrict our 
attention to the case of coherent scattering in the laboratory frame. The QED theory 
of dipole scattering, for the more realistic case of complete frequency redistribution 
(CRD) and partial frequency redistribution (PRD) are respectively developed by 
Landi Degl'Innocenti & Landolfi ([28], and references cited therein) and by Bommier 
[35[ [36l [37[ [38] . These authors consider Hanle scattering in the presence of only 
magnetic fields of arbitrary strength. The classical non-perturbative theory of Hanle- 
Zeeman scattering with PRD in atom's rest frame is presented in [39]. Sampoorna et 
al. [^ Hi] present the corresponding laboratory frame expressions in a form suitable 
for astrophysical applications. 

A classical theory for the coronal forbidden emission lines that arise from magnetic 
dipole (Ml) transition is derived in [9] and [lO]. The purpose of the present paper 
is mainly to derive the phase matrices characterizing atomic scattering of polarized 
radiation for forbidden transitions in external magnetic fields. This approach can handle 
all multipoles in the transition. In particular, we present the explicit forms of the phase 
matrices for E2 and M2 type transitions, as they are not available in the literature. In 
sections 2 and 3 we present the theoretical framework briefly. In section 4 we describe 3 
limiting cases of scattering in magnetic fields. The numerical examples illustrating the 
application of these formulae are given in section 5. 

2. Theory of scattering for forbidden transitions 

It is well-known that an energy level of an atom with total angular momentum J 
splits into (2J + 1) equally spaced components (Zeeman splitting) when the atom is 
exposed to an external magnetic field B. These levels are characterized by states | Jm), 
where m denotes the magnetic quantum number which is projection of the total angular 
momentum operator J along B, which is taken along the Z— axis. If we consider atomic 
scattering Jj — J — >■ of radiation, the elements of the on-energy-shell transition 
matrix T are given [311 132] by 

{Jfmf\T\Jimi) = y^^{Jfmf\S{k, fi)\Jm)(prn{Jm\A{k' , fi')\Jimi) , (1) 
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where (fm denotes the profile function 

V^m = i^mnif — UJ — iT) ^ ] ^mnif = — Ef . (2) 

Here u = 2t{v denotes the frequency of the scattered radiation, Ei, Em and Ef denote 
respectively the energies associated with the initial |Jjmj), intermediate \Jm) and 
final \JfiTLf) states of the atom, and F denotes the natural width associated with the 
intermediate state. The vectors k' {k', 9', (f)') and k {k, 9, (p) are photon momenta of the 
incident and scattered radiation. The states of circular polarization associated with k' 
and k are denoted respectively by = ±1 following the convention employed by 
Rose [l2]. If p{k') denotes the 2x2 density matrix representing the state of polarization 
of the incident radiation, the density matrix of the scattered radiation is given by 

p{k) = Tp{k')T\ (3) 

where T''^ denotes the hermitian conjugate of T considered as a 2 x 2 matrix with respect 
to the basis states of circular polarization. Following [311 [32l [331 [M] , the matrix elements 
for emission and absorption processes in Eq. ([l]) are given by 

{Jfmf\£(k,ii)\Jm) = ^C(JyLJ;mjMm)D^,^(0,^,O)*(-i/i)'^(^)jiH*, (4) 

L 

{Jm\A{k\^l')\Jm^) = J2C{J,L'J; m,M'm)D^j,^,{<P',9',0){-ifiT'^'''^jL'{uj'), (5) 

L' 

where uj = k, uj' = k' and denotes the reduced matrix elements for the atomic 

transition from a lower level to an upper level given by Eq. (31) of [31]. Further 

h{L) = i[l + 7r;7r(-l)^] ; h{L') = ^[1 + 7r,7r(-l)^'], (6) 

with TTj, TT and vr/ denoting respectively the parities of the initial, intermediate and final 
atomic states. All transitions other than L = L' = Tnif = UTTi = 1 are referred to as 
forbidden, following the traditional usage. 



3. Hanle-Zeeman scattering phase matrix 

Any density matrix p for polarized radiation may be expressed in terms of the 
corresponding Stokes parameters Sp with p = 0,1, 2, 3, through 

1 ^ 

P=2 5Z^A. (7) 

where (Tq denotes the 2x2 unit matrix and <Ti^2,3 denote the well-known Pauli matrices. 
The parameters {Sq, Si, S2, S3) = {I ,Q,U,V) denote the Stokes vector S. The 4x4 
scattering matrix 72. is now defined through the relation 

S{k)=nS'{k') . (8) 

Noting that 

Sp = tr{app), (9) 
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where tr denotes the trace, the elements of the scattering phase matrix are given by 

T^PP' = ltr{a,Tcr,, Tt) = 1 ' , (10) 

mm' 

where the 4x4 matrix "P™™ is the phase matrix for hne scattering, which is also called 
'scattering matrix' in the literature. The product <^m fm' can be converted into a sum 
using the conversion formula of [l3[ |39] , to which a Doppler convolution is applied. In 
scattering on a two-level atom (Jj = Jf = 0), the elements of ■p™'" can be expressed in 
an elegant form 

P-r' = g;'"''{ky g^^'^'ik') , (ii) 

where 

X J2 CiLLl; m - m'mi)C{LLl- - /iVz)/^^,^, (0, 0, 0) , (12) 

and {k') is given by a similar expression where p, L, u are now replaced by p', L', uj' 
respectively. Equation (fTOl) can be used to consider scattering in magnetic fields of ar- 
bitrary strength and orientation. Therefore we refer to 72. as the Hanle-Zeeman phase 
matrix. It is interesting to consider the particular case of ^ 2 — scattering for 
forbidden transitions with L = L' = 2 and study the phase matrix. We consider three 
domains of the field strength B, namely 

(i) the case where B is sufficiently strong so that the levels with different m are distinct, 
leading to "Zeeman scattering" with m = m' as the summation over m in Eq. ffTOl) drops 
outs (no m-state interference) ; 

(ii) the weak field limit, where the magnetic splitting of the upper level is of the same 
order as the natural width F, so that the quantum interferences between magnetic sub- 
levels m and m' take place leading to the well-known "Hanle scattering" ; and 

(iii) "Rayleigh scattering" in the absence of magnetic field, when the levels with different 
m are degenerate, so that there is only a single phase matrix 'P'"™' = "P. 

In the Hanle scattering case, one has to consider totally 25 pairs of m and m' for 
each element of the 4x4 phase matrix p™-™- given by Eqs. (iTTl [T2l) with L = 2 and 
m = m' = 2,1,0,-1,-2, while for the Zeeman scattering, one has 5 phase matrices 
pmm ^ fpm ^j^-]^ ^ _ _|_2^ -|-]^^ ^^^e prcscut thcsc 3 limiting cases in the following 
sections. 
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4. Limiting Cases 



4.I. Zeeman scattering phase matrix (Rayleigh scattering in strong fields) 



The Zeeman scattering matrix for J = 
can be written as 

5 



n 



Zeeman 



transition in the strong field hmit 



(13) 



where V"^ and H{vm, a) represent the scattering phase matrices and the Voigt profile 
functions containing the energies of the upper states \Jm) with m = —2, ■■■,2 
respectively. The normalization constant is 5/4. The frequency dependent shift Vm 
and the damping parameter a in Voigt profile function are defined following |43j . 

Zeeman scattering phase matrices V"^ for electric quadrupole {E2) and the magnetic 
quadrupole (M2) forbidden line transitions are given below. In the matrix elements, the 
upper sign represents M2 transition and the lower sign an E2 transition. M2 and E2 
transition have the same expression for some matrix elements (eg. the (1, 1) element). 



p2 = ±(l-5'2)(l_52) 



/ |(l + <5'2)(i + 52) ±l{i-5'^)(i + S^) <5'{l + 52) \ 




V <5(l + <5'2) 



±5(1 -5'2) 



255' 



(14) 



-(l-(5'2)(l-52^ 



f i(l+<5'2)(l+52) ±l(l_y2)(i^52) _y(l+52)^ 

±^(l + <5'2)(l-52) l(l_5'2)(i_52) ^5'(l-52) 



V -(5(l + <5'2) ^5(i_5'2) 255' J 



(15) 



/ ^[(l-25'2)2 +5'2][(i_ 252)2^^2] ^ 1 [(1 _ 2^2)2 _y2][(i_ 2^2)2.^ ^2] q qs^ 

T^[(l-25'2)2+5'2][(i_ 252)2 -52] l[(l_25'2)2_y2][(^_ 2^2)2 _52] 



oj 



V -i5(l -252)[(i_25'2)2 + 5'2] 



±i5(l-252)[(i-25'2)2_5'2] 



/O -i5'(l-25'2)[(l -252)2 + 52] 
±i5'(l-25'2)[(l- 252)2 - 52] 





Vo 



(16) 



55'(l-25'2)(l- 252) / 



/' i[(l-25'2)2+ 5'2][(1 - 252)2 + 52] +i[(i-25'2)2-5'2][(l- 252)2 + 52] 

T^[(l-25'2)2+5'2][(l- 252)2 - 52] 1[(^_25'2)2-5'2][(1 - 252)2 - 52] 



V i5(l - 252)[(1 - 25'2)2 +5'2] +i5(l- 252)[(l-25'2)2-5'2] 0/ 



/O |5'(1 -25'2)[(1 - 252)2 + 52] 

T35'(l - 25'2)[(1 - 252)2 - 52] 



Vo 55'(1 -25'2)(1 - 252) J 



(17) 



V° = 95'^5'^{1-5'^){1~5'^) 



/I Tl 0\ 

=Fl 



V 0/ 



(18) 
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where the symbols 5' = cos 6' and 5 = cos 6, with 6' and 6 denoting the angle made 
by the incoming and outgoing ray with respect to the magnetic field chosen to be along 
the Z-axis of a coordinate system. 

If we average over all the incoming angles, we obtain the emission coefficients f]i,Qy 
used in the thermal radiation limit (LTE). rju = for the particular choice of the 
coordinate system used in this calculation. However, rju can be generated by a clockwise 
rotation of the coordinate system through an azimuthal angle x- We have verified that 
the new emission coefficients rjj^Q^uy obtained after rotation are identical with those 
derived by Beckers [13] for the quadrupole emission in Zeeman effect regime. 



4-2. Rayleigh scattering phase matrix (the non-magnetic case) 

The Rayleigh scattering matrix for J = ^ 2 ^ transition in the absence of magnetic 
fields is given by 



n 



Rayleigh 



-VH{v,a) 



(19) 



In this case all the magnetic sublevels of the J = 2 are degenerate. Thus 
the profile functions (pm with m = —2,-1,0,1,2 become identical and equal to ipo. 
Therefore we obtain only one scattering phase matrix V, and the profile function H{v, a) 
corresponding to ipo. The elements of the Rayleigh scattering phase matrix are given by 



Poo 



-(1 - 35^ + i5*){l - 3<5'2 + 45"') + (1 - 5^){1 - 5'^) 



+ ci(l - 52)1/2(1 _ 5'2)i/2 [4535/3 _^ g5y(^ _ 252)(i _ 25'2)] 
+ C2 6(1 - 5'')(1 - 5'2)52 5'2 + i(l - + 45'*)(1 - 55'^ + 45'*) 

+ 4C3(1 - 52)3/2(1 _ y2)3/25y ^ 1^^(^ _ g2-^2^-^ _ y2)2 ^ 



(20) 



■Poi 



::p _(1 - 352 + 45*)(1 - 55'2 + 45"*) + (1 - 52)(1 - 5'2) 



±i(l + 52)(i_ 5/2-) ^9525/2 



+ ci(l - 52)1/2(1 _ 5'2)i/2 [±4535'(1 - 5'2) =F 655'(1 - 252)(1 - 25'2)] 



T C2 



(1 - 52)(1 - 5'2)52 5'2 + 1(1 - 552 + 454)(1 - 35'2 + 45'*) 



± 4C3(1 - 52)3/2(1 _ 5'2)l/25y3 ^ 1 (^ _ ^2^2^^ _ g,2~^^^ ^ y2) 



(21) 



Po2 = ±si(l- 52)1/2(1 -5'2)i/2 [253(1 -5'2) + 655'2(l- 252)] 
+ S2 [±5'(1 - 25')(1 - 552 + 45*) T 6(1 - 52)(1 - 5'2)52 5'] 
=F 2S3(1 - 52)3/2(1 _ 5'2)1/25(1 - 35'2) ± S4(l - 52)2(1 - 5'2)5' 



(22) 



Vio= T ^(1 - 352 + 45*)(1 - 55'2 + 45'*) + (1 - 52)(1 - 5'2) 

■5'3(1 - 52) =F 655'(1 - 252)(1 - 25'2)] 

I 



+ ci(l- 52)1/2(1 -5'2)i/2 [±455 

1 - 5'2)52 5'2 + 1(1 - 352 + 45*)(1 - 55'2 + 45'*) 



l(l-52)(l + 5'2):f 9525'2 



T C2 



6{l-52)(: 

± 4C3(1 - 52)1/2(1 _ 5'2)3/253y ^ l^^(i _ 52)(^ _ y2)2(^ ^ ^2) 



(23) 
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Vii 



1 



(1 - 55^ + 45'^){1 - 5<5'2 + 45"') + (1 - 5^){1 - S''^) 



{1 -5^) {1-5'^) + 95^5'- 



+ ci(l - 52)1/2(1 _ 5'2)l/2 [45^(1 _ ^2-)(i _ _^ g^y^i _ 252)(1 _ 2^2)] 



+ C2 



6(1 - <52)(1 - 5'2)52 5/2 _^ _ ^ 454)(i „ 3^2 4^4) 



+ 4C3(1 - 52)1/2(1 _ 5'2)l/2^3y3 _^ i ^_^(i _ 52)(i _ y2)(i ^ 52)(i _^ y2) ^ 



(24) 



Pi2 = si{l - 52)1/2(1 _ y2^i/2 [25(1 _ s^){l - 5'2) - 655'2(1 - 252)] 

- S2 [(1 - 352 + 454)5'(1 - 25'2) - 6(1 - 52)(1 - 5'2)52 5'] 

- 2s3(l - 52)1/2(1 - 5'2)i/253(l - 35'2) + S4(l - 52)(1 - 5'2)5'(1 + 52) , (25) 



V20 = T ^1(1 - 52)1/2(1 _ 5'2)l/2 [2(1 _ 52)^3 _^ g^2y(i _ 25'2)] 
+ S2 [±6(1 - 52)(1 - 5'2)55'2 =F 5(1 - 252)(1 - 55'2 + 45'^)] 

± 2S3(1 - 52)1/2(1 _ 5'2)3/25'(l _ 3^2) ^ 5^(1 _ 52)(i _ ^/2)2^ ^ (36) 



V21 = - si{l - 52)1/2(1 _ 5'2)i/2 [25'(1 - 5'2)(1 - 52) - 6525'(1 - 25'2)] 
+ S2 [5(1 - 252)(1 - 35'2 + 45"i) - 6(1 - 52)(1 - 5'2)55'2] 

+ 2S3(1 - 52)1/2(1 - 5'2)1/25'3(1 - 352) - S4(l - 52)(1 - 5'2)5(1 + 5'2) , (27) 



V22 = ci(l - 52)1/2(1 - 5'2)i/2 [(1 - 52)(1 - 5'2) + 652 5'2] 
+ C2 [6(1 - 52)(1 - 5'2)55' + 255'(1 - 252)(1 - 25'2)] 

+ C3(l - 52)1/2(1 _ 5'2)i/2(l _ 352)(1 - 35'2) + 2c4(l - 52)(1 - 5'2)55' , (28) 



P33 = 255'(1 - 52)(1 - 5'2) + 255'(1 - 252)(1 - 25'2) 

+ ci(l - 52)1/2(1 _ 5'2)i/2 [(1 - 352)(1 - 35'2) + 6525'2] 

+ 6c2(l - 52)(1 - 5'2)55' + C3(l - 52)1/2(1 - 5'2)i/2(l - 52)(1 - 5'2) , (29) 
7^03 = 7^30 = Vvj. = •P31 = P23 = P32 = . (30) 

In the above equations, the symbols c„ = cos?2(0 — 0') and s„ = sinn(0 — 0'). The 
incoming and outgoing rays make an angle {6', 0') and {6, 0) with respect to the Z-axis 
of a coordinate system. 

4-3. Hanle scattering phase matrix (in weak field limit) 

In the weak-field limit, the magnetic splitting is of the same order as the natural line 
width of an atomic state. Therefore there exists a frequency coupling between the 
magnetic sublevels of the upper energy level. If the product of profile functions (fm^m' 
is simplified as mentioned before (see the paragraph following Eq. fllUI) ). the weak field 
Hanle phase matrix is obtained. 
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5. Numerical results for a — > 2 ^ transition 

We consider the example of single scattering involving a — > 2 — > forbidden transition. 
An unpolarized beam of radiation is incident on the atom. Therefore the four Stokes 
parameters {I,Q,U,V) of the scattered radiation are (7?.oo, "^lo, "^20, "^so) respectively. 
In the expressions ( 1311) - ( 137|) for the Stokes parameters, the upper sign represents M2 
type transition and the lower sign represents E2 type transition. The different curves in 
all the figures represent difi^erent values of magnetic field strength defined by a splitting 
parameter vb = (^z/l/Az/d, where z/l is the Larmor frequency, g the Lande factor, and 
Avd the Doppler width. The solid curve corresponds to vb = 0.004, the dotted curve 
to Vb = 0.02, the dot-dashed curve to 0.1, the dashed curve to 0.5 and the long-dashed 
curve to 2.5. The damping parameter of the Voigt profile function a = 0.004. The 
direction of the magnetic field is chosen along the Z-axis of the co-ordinate system. 

5.1. Forward scattering 

Figure [1] shows the scattered Stokes profiles in a forward scattering event {6' = 0°, 0' = 
0°; 6 = 0°, (f) = 0°). For this choice of scattering geometry, the Stokes parameters, given 
by the following expressions, are the same for both M2 and E2 type transitions. For 
the case of non-magnetic Rayleigh scattering : 

/=^boV^o]; Q = 0; f/ = 0; V = 0. (31) 
Likewise, the scattered Stokes profiles for the Hanle-Zeeman scattering are : 

i = l[^i^l + ^-i^U]; Q = 0; u = o- v = ^[^,^i-y,^,^u]- m 

In the particular geometry of forward scattering, the same expressions hold good 
for the Zeeman scattering case also. Therefore we show only / and V/I profiles 
corresponding to the Hanle-Zeeman case. The Stokes I for the Rayleigh scattering 
case is nearly same as the Stokes I of the Hanle-Zeeman case (the solid curve). The 
Stokes Q/I = in both the cases. It can be understood from the classical theory that 
the degree of linear polarization is zero for the forward Rayleigh scattering. As the 
line of sight is taken parallel to the magnetic field, the standard doublet pattern of 
longitudinal Zeeman effect is seen in /. The V/I profiles show anti-symmetric pattern 
(oppositely circularly polarized). These components are produced due to the scattering 
involving |0,0) and |2,±1) states. The resolution of / into 2 components takes place 
only for very strong fields like vb = 2.5. For weaker fields, it appears only as a slight 
broadening of the non-magnetic / profiles. All the weak field cases for vb < 0.1 merge, 
and are not distinguishable graphically. 

5.2. Scattering at 90° 

This scattering geometry corresponds to the maximum degree of linear polarization. 
The angles chosen for this case are : 6' = 0° , (p' = 0° and 6 = 90°, = 45°. In this case. 
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the scattered Stokes profiles are given by 

/ = ^bo^o]; Q = T^bo^o]; ^ = 0; ^ = 0, (33) 

for Rayleigli scattering case, and 

I = l[^^^l + v9*_i] ; Q = tI [Vi + V^i V-i] ; U = 0; V = 0, (34) 

o o 

for the Hanle-Zeeman as well as the Zeeman scattering cases. 

It is clearly seen that the maximum degree of linear polarization Q// = —1 in all 
the three cases. The Q/I ratio is independent of frequency as well as the magnetic field 
strength. From Eqs. fl33p and (1341) . we observe that the linear polarization (Q) profiles 
for M2 and E2 type transitions differ only in sign. As the I profiles look very similar 
to the I profiles of Fig. [T], except for the fact that U/I = V/I = and Q/I = —1, we 
do not show these profiles again. 

5.3. Scattering at an arbitrary angle 

In Figs. [2] and [3] we show the results for M2 and E2 type transitions respectively. They 
correspond to the choice of scattering geometry : 6' = 45°, (f)' = 0° and 6 = 90°, (p = 45°. 
For this specific choice of angles, the scattered Stokes parameters are given by the 
following simple expressions : 

1 = 1 [fofo] ; Q = 0; U = 0; V = 0, (35) 
for non-magnetic Rayleigh scattering case. For the Zeeman scattering case we have 

I = ^ [0.125 ((^i(^t + (^_iV9*_i) + 0.1875 ((^2<^; + (^_2<2)] , 

Q = ^ [tO.125 (ifM + if-iifU) ± 0.1875 (y^a^^ + ip-2ipU)\ , (36) 

with U = = 0, and finally for the Hanle-Zeeman case we get 
5 

/ = - [0.125 {^i^l + i^-ivU) + 0.1875 {^2^*^ + 

+ i X 0.125 {^-i^l - VivU) - 0.0625 {^2^*_2 + ^-2^*2% 
g = ^ [ T 0.125 {^i^l + V-ivU) ± 0.1875 {^2^1 + ^-2VU) 

T i X 0.125 {^-i^l - Vi^U) T 0.0625 {^2V*_2 + ^-2^1)1 
[/ = ^ X 0.0883883 [(1 - i) {±^-2^1 T ^2^1 ± V-iVl T ^-ivU) 

+ (1 + i) {±ipnp*_2 =F V1V2 ± V2V*^i T V-2V-i)]^ 
V 0.0883883 [(1 - i) - ^i^l + V-2^-1 - ^2^-1) 

+ (1 + i) {ip-2V*i - ^2^1 + V-1^-2 - ^-1^2)]- (37) 

Clearly the Q/I and U/I in Figs. [2] and [3] for both Zeeman as well as Hanle-Zeeman 
cases differ only in sign (see Eqs. fl36p and fl37p ). However the V/I remains the same 
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for both M2 and E2 type transitions. In the following we discuss only Fig. [2] for M2 
type transition. 

In the cases of Zeeman and Hanle-Zeeman scattering, the Stokes / profiles clearly 
show splitting in the strong fields. When the fields are weak, these components are 
not resolved. In this choice of geometry, the Stokes Q gets contributions not only from 
the energy states |2,±1) but also from |2,±2) in both the cases (see Eqs. (|36|) and 
fl37|) ). Rayleigh scattered Stokes Q is zero in the M2 as well as E2 transition, even 
for a non-zero scattering angle (see Eq. ( 135|) ). The / profile for Rayleigh scattering is 
similar to the solid curve of the general Hanle-Zeeman case, and hence is not shown 
again. Equations (l36l) and (1371) . differ only in the m-state interference terms. The effect 
of these extra terms in Eq. fl37j) becomes significant only for weak fields (compare solid, 
dotted and dot-dashed Q/I curves in the 2nd row of Fig. [2l) . When the fields are stronger 
(dashed and long-dashed curves), the Zeeman and Hanle-Zeeman scattering give nearly 
same Q/I profiles, showing that the interference terms gradually become negligible. 

The Stokes U and V are zero in Rayleigh and Zeeman scattering cases due to our 
choice of the scattering geometry, but they appear in Hanle-Zeeman scattering cases in 
the weak field limit. The generation of these U and V in the Hanle-Zeeman scattering 
case is entirely due to the interference between different m-states. As the magnetic field 
strength increases, U/I decreases and approaches small values. The V/I is extremely 
small in the weak field limit and gradually develops the typical anti-symmetric profiles 
of the usual Zeeman effect as the field strength increases. 

6. Conclusions 

The polarization phase matrices for forbidden line transitions are derived. The 3 
important regimes of Rayleigh, Zeeman, and Hanle-Zeeman scattering are studied. The 
regime of Hanle-Zeeman scattering provides a complete description of scattering theory 
in magnetic fields of arbitrary strengths, with the other two regimes serving as extreme 
limiting cases. We show that the linear polarization profiles {Q and U) for M2 and E2 
type transitions differ only in sign. The V profiles are however, insensitive to the type 
of transition. 
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Figure 1. The Stokes profiles in the regime of Hanle-Zeeman scattering. The case of 
forward scattering is considered. Different curves represent various values of the field 
strength vb- See sections [5l and [Ql for model parameters. 
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Figure 2. Zeeman (l*** column) and Hanle-Zeeman scattering (2"*^ column) in a 
M2 type scattering event. The scattering geometry is defined through the choice of 
angles : B' — 45°, = 0° and G ~ 90°, — 45°. Different curves correspond to different 
values of the field strengths expressed through the parameter vb- The solid curves: 
vb = 0.004, dotted curves: 0.02, dot-dashed curves: 0.1, dashed curves: 0.5, and long 
dashed curves: 2.5. See section [Q] for discussions. 
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Figure 3. Zeeman (I*'* column) and Hanle-Zeeman scattering (2"^^ column) in a E2 
type single scattering event. The scattering geometry and different curve types are the 
same as in Fig. [2j 



